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Abstract 

Uniaxial  tension  tests  using  a  strain  rate  of  0.04  in./in./min.  are  performed  on  rectangular  smooth 
and  single  edge-notched  specimens  of  varying  thicknesses  for  a  composite  solid  propellant. 
Stress-strain,  crack  growth,  crack  growth  rate  and  crack  growth  resistance  data  are  provided. 
Thickness  effects  and  the  mechanism  of  crack  growth  are  described.  Methods  of  calculation  are 
explained  for  the  crack  growth  rate  and  the  Mode  1  stress  intensity  factor.  A  model  is  developed 
for  the  crack  stable  growth  rate  as  a  function  of  the  stress  intensity  factor. 

Key  Words 

linear  elastic  fracture  mechanics,  LEFM,  plane  stress,  composite  solid  propellant,  crack  growth, 
crack  growth  rate,  mode  1  stress  intensity  factor,  crack  growth  resistance  curve,  thickness  effect. 

Introduction 

Composite  solid  propellants  consist  of  a  lightly  crosslinked  polymer,  highly  filled  with  relatively 
coarse  solid  particles.  They  are  considered  to  be  incompressible.  Their  loaded  response  is 
viscoelastic.  Previous  studies  have  been  conducted  to  investigate  crack  growth  in  these  materials 
{1-5}.  The  basic  approach  is  based  on  linear  elastic  or  linear  viscoelastic  fracture  mechanics. 
Experimental  results  have  shown  that  a  power  law  relationship  exists  between  the  crack  growth 
rate  (a)  and  the  Mode  1  (opening  mode)  stress  intensity  factor  (K^ .  This  is  consistent  with  the 
linear  viscoelastic  fracture  theories  developed  by  Knauss  {6}  and  Schapery  {7}. 

An  example  of  a  composite  solid  propellant  can  be  described  as  follows  {8}.  Ammonium 

perchlorate  crystal  particles  of  two  sizes,  diameters  20  p  (lp  =  10-6  meter)  and  200p  function 
as  the  oxidizer.  The  oxidizer  weight  fraction  for  the  20p  particles  is  0.3.  The  oxidizer  weight 
fraction  for  the  200p  particles  is  0.7.  Aluminum  powder  particles  of  diameter  30p  function  as 
the  fuel.  Polybutadiene  rubber  functions  as  the  matrix/binder  (also  as  fuel).  The  weight  fractions 
for  the  oxidizer,  fuel,  and  matrix/binder  are  respectively  0.7,  0.2,  0.1.  The  specific  gravities  are 
respectively  1.69,  2.65,  1.0.  The  volume  fractions  are  therefore  respectively  0.7,  0.13,  0.17.  The 
particle  volume  fraction  (PVF)  is  0.83.  For  every  200p  diameter  particle,  there  are  approximately 
77  30p  diameter  particles,  and  429  20p  particles. 

For  this  paper,  four  different  thicknesses  (0.2  in.,  0.5  in.,  1.0  in.,  1.5  in.)  of  smooth  and  single 
edge-notched  tension  (SENT)  specimens  of  a  composite  solid  propellant  were  uniaxially  strained 
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at  a  rate  (e)  of  0.04  in./in./min.  The  induced  end  force  (F)  and  crack  length  (a)  (for  the  SENT 

specimens)  were  recorded  as  functions  of  e  (equivalently  as  functions  of  time,  t ).  Photographs 
showing  the  crack  growth  mechanism  are  included  in  the  paper,  and  the  mechanism  is  described. 
Specimen  stress  (a)  is  plotted  versus  specimen  strain  (e)  for  both  the  smooth  and  notched 
specimens,  a  is  calculated,  and  a  and  a  are  plotted  versus  t.  Kx  is  calculated  and  plotted  versus 

the  change  in  the  crack  length  (8a) ,  giving  the  material’s  crack  growth  resistance  curve  (R- 
curve).  There  are  plots  showing  the  relationship  between  a  and  Kx .  The  a(Kx)  function  is 

modeled  and  interpreted.  Effects  of  specimen  thickness  (b)  are  noted.  Linear  elastic  fracture 
mechanics  (LEFM)  was  used.  For  finite  element  calculations  the  state  of  plane  stress  was  used  for 
all  thicknesses.  This  is  because  there  was  no  consistent  evidence  of  plane  strain,  even  in  the 
thickest  (1.5  in.)  specimens,  as  will  be  reported.  Because  of  the  material’s  assumed 
incompressibility,  0.499  was  used  as  the  value  of  Poisson’s  ratio  (v). 

Experiments 

All  test  specimens  were  rectangular,  1  in.  wide  (w  =  1  in.,  perpendicular  to  the  direction  of 
loading)  and  5  in.  long  (h  =  5  in. ).  Specimen  ends  were  glued  to  aluminum  tabs  and  subjected  to 
an  applied  uniform  uniaxial  displacement  rate  (A)  of  0.2  in./min.  (applied  uniform  uniaxial 

nominal  strain  rate  (e)  of  0.04  in./in./min.).  Four  specimen  thicknesses  (b)  were  used,  0.2  in.,  0.5 
in.,  1.0  in.,  and  1.5  in.  Eight  unnotched  specimens  were  tested,  two  of  each  thickness.  The  stress 
(a)  vs.  strain  (e)  curves  for  those  tests  are  shown  in  Figure  7.  The  average  initial  (linear)  value 
(E0 )  for  Young’s  Modulus  is  734  psi.  The  average  maximum  average  specimen  stress  (cmax  )  is 
90.9  psi,  and  the  average  strain  at  the  maximum  stress  (£(amax) )  is  0.187.  The  average  maximum 
strain  (emax  )  is  0.204.  The  overall  trend  is  for  (emax )  to  decrease  as  b  increases.  Eight  single 
edge-notched  (SENT)  specimens  were  tested  (Mode  1  loading),  two  of  each  thickness.  The  notch 
length  (aQ ,  cut  along  the  specimen  width)  used  was  0.3  in.  The  o  vs.  e  curves  for  those  tests  are 

shown  in  Figure  8.  (a  is  calculated  using  the  specimen  uncracked  cross  sectional  area,  A  =  bw). 
The  curves  show  small  slopes  near  the  origin.  This  is  probably  because  the  specimens  were 
unintentionally  put  in  an  initial  state  of  slight  compression.  The  average  value  for  amax 

(calculated  using  A  =  bw)  is  40.6  psi.  The  average  value  for  £(crmax)  is  0.075.  The  overall  trend 
is  for  Young’s  Modulus  ( E )  and  Gmax  to  increase  as  b  increases. 

The  plots  of  Figures  8-14  refer  to  the  eight  SENT  specimens.  In  each  plot  the  data  shown  for  a 
particular  specimen  stops  when  failure  begins,  and  therefore  corresponds  to  the  strain  range 
0  <  £  <  £(amax)  •  This  is  done  so  that  LEFM  can  be  used.  LEFM  is  not  valid  once  a  begins  to 

decrease  as  a  specimen  fails. 

Crack  Growth  Mechanism 

The  largest  particles  in  a  composite  solid  propellant  are  typically  larger  than  100(X  (0.1  mm.),  so 
at  that  scale  these  materials  are  non-homogeneous.  When  stretched,  the  variations  in  particle  size, 
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particle  distribution,  cross  link  density,  and  particle-binder  bonds  result  in  non-homogeneous 
stresses  and  strengths.  Because  the  particles  can  be  considered  rigid  with  respect  to  the  polymer, 
the  magnitudes  of  local  stresses  can  be  much  greater  than  the  magnitudes  of  the  applied  stresses, 
particularly  when  the  gaps  between  particles  are  small.  Because  of  the  randomness  of  local 
stresses  and  strengths,  local  failure  locations  usually  do  not  coincide  with  the  points  of  maximum 
stress  as  determined  from  an  elastic  analysis  for  a  homogeneous  isotropic  material.  For  a  cracked 
specimen,  a  damage  zone  forms  ahead  of  the  crack  tip.  Local  failures  (voids  or  microcracks) 
usually  occur  first  in  the  interior  of  the  damage  zone,  not  at  the  crack  tip  itself.  The  crack  tip 
extends  into  the  damage  zone  when  the  ligament  of  material  separating  it  from  the  nearest  void  or 
microcrack  breaks.  The  front  of  the  void  or  microcrack  becomes  the  new  crack  tip.  As  the  crack 
tip  extends,  the  damage  zone  moves  forward,  gaining  new  material.  The  type  of  damage  which 
causes  void  formation  is  matrix/binder  cracking.  The  types  of  damage  which  cause  damage  zone 
microcracking  are  particle-matrix/binder  debonding  (dewetting)  and  particle  cracking. 

Figures  1-6  show  photographs  of  crack  growth  in  the  SENT  specimens.  All  four  specimen 
thicknesses  (0.2  in.,  0.5  in.,  1 .0  in.,  1 .5  in.)  are  shown,  and  the  common  growth  mechanism  is  void 
formation  in  the  damage  zone  ahead  of  the  crack  tip,  and  then  crack  tip  extension  due  to  ligament 
rupture.  This  type  of  growth  produces  the  rough  crack  surfaces  shown  in  Figures  3-6.  Figures  1 
and  2  show  crack  initiation  from  an  initial  notch.  Crack  tip  blunting  occurs  before  initiation 
because  of  the  large  uniaxial  strain  (Figure  7)  the  material  is  capable  of  withstanding. 

Examinations  of  the  SENT  specimen  fracture  surfaces  showed,  that  for  all  thicknesses  the  crack 
fronts  were  straight  lines  and  perpendicular  to  the  direction  of  crack  growth.  This  implies  that  the 
region  near  the  crack  tip  always  experienced  a  state  of  plane  stress.  If  plane  strain  existed,  it 
would  have  been  indicated  in  the  thicker  specimens  by  a  curved  crack  front  which  was  advanced 
farther  near  its  middle,  and  less  near  its  edges.  In  that  case  the  interior  of  the  crack  tip  region 
would  be  in  plane  strain  and  the  exterior  would  be  in  plane  stress.  (A  material’s  crack  growth 
resistance  is  less  in  plane  strain  than  in  plane  stress).  Apparently,  void  formations  in  the  damage 
zone  prevented  the  lateral  constraint  which  is  necessary  for  a  state  of  plane  strain. 

Calculation  of  the  Crack  Growth  Rate 

The  rates  (a)  of  crack  growth  (for  n  +  1  points  (a,  t )  numbered  0, 1, ...,  n)  were  calculated  by 


m  = 


a(l)  -c(0) 
f(l)-r(0) 


a(i)  =  max 


a(j)  -  aji  -  1)  a(i  +  1 )  -a(Q') 
1)’  t(i+  l)-t(i)  ) 


a(n )  = 


ajn)  -a(n  -  1) 
t(n)-t(n-  1) 


(1) 


This  method  is  easy  to  implement,  and  assigns  each  experimental  point  (a,  t)  a  value  of  a .  The 
values  of  a  are  not  smoothed.  Each  interior  point  (a,  t )  is  assigned  the  larger  of  the  two  adjacent 
calculated  values.  Using  the  method  results  in  the  ability  to  determine  the  most  conservative 
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upper  bound  on  the  envelope  of  a  for  the  set  of  SENT  specimens. 

The  Experimental  Crack  Growth  Rate 

Figure  9  shows  a  plot  of  crack  length  (a)  vs.  time  ( t ) .  The  crack  initiation  time  (tc)  varies  from 
1.34  min.  to  1.59  min.  The  average  value  for  tc  is  1.45  min.  The  corresponding  applied  nominal 
strain  at  crack  initiation  (ec  =  tc  •  e)  varies  from  0.0535  to  0.0636.  The  average  value  for  ec  is 

0.0581.  The  maximum  recorded  value  for  a  is  0.427  in.  Figure  10  shows  a  plot  of  crack  growth 
rate  (a)  vs.  time.  The  maximum  calculated  value  for  a  is  0.847  in./min.,  approximately  4.2  times 
the  applied  displacement  rate  of  A  =  0.2  in./min. 

Calculation  of  the  Mode  1  Stress  Intensity  Factor 

For  the  SENT  experiments,  a  specimen’s  edge  displacement  (A)  was  prescribed  as  a  function  of 
time  (r).  The  crack  length  (a),  was  measured  as  a  function  of  t.  For  each  point  (A,  a)  of  each 
specimen  thickness  ( b ) ,  the  LEFM  Mode  1  stress  intensity  factor  (  K  ^  )  was  calculated  by  the 
method  described  in  this  section. 

K  j  can  be  written  in  the  form 


*1 


(2) 


where  c£vg'  is  the  average  value  of  the  specimen  edge  stress  (oM)  induced  by  the  application  of 
the  uniform  edge  displacement  (A) ;  w  is  the  specimen  width,  h  is  the  specimen  length,  and  /  is 
some  function  of  a/h  and  h/w .  Dividing  by  o™g'Ja  gives 


K, 


O*  avg.  ^ 


(3) 


where  g  is  a  second  function  of  a/h  and  h/w.  According  to  LEFM,  the  energy  release  rate  (J) 
(for  plane  stress)  can  be  expressed  as 


J  = 


(4) 


where  E  is  Young’s  Modulus.  (For  the  case  of  plane  strain,  in  Equation  4,  E  —>  E/(  1  -  v2) .)  This 
implies  that 
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K{  =  4je. 


(5) 


olvg'  can  be  expressed  as 


(6) 


where  F  is  the  end  force  induced  on  the  specimen,  and  A  is  the  specimen’s  uncracked  cross- 
sectional  area  (bw) .  Therefore  g'  becomes 


Using  Equations  5  and  7  in  Equation  3  gives 


(7) 


(8) 


The  function  g( —  ]  is  the  specimen  geometry  correction  factor  for  stress  intensity,  and  is 
\w  wj 

1  f  J  TJI 

approximated  by  numerically  calculating  -pr  J—  *  Linear  elastic  plane  stress  finite  element 

calculations  (ABAQUS  code)  were  performed  of  a  SENT  specimen  subjected  to  a  uniform 
uniaxial  displacement  (A)  ,  using  w  =  1  and  h  =  5  .  The  value  of  the  specimen  thickness  ( b ) 
does  not  appear  in  plane  finite  element  calculations,  and  is  implicitly  equal  to  1 .  (Note  that  the 
value  of  g  is  independent  of  the  value  of  the  prescribed  displacement  (A)  used  in  a  finite  element 

calculation.  Note  also  that  although  E  appears  in  the  expression  f°r  8  > the  value  °f the 

expression  is  independent  of  the  value  of  E  used.  Any  (positive)  values  can  be  used  for  A  and  E 
in  the  finite  element  calculations.  The  value  of  g  depends  only  on  the  two  parameters  a/ w  and 
h/w .)  The  value  of  a  was  varied,  and  one  finite  element  calculation  was  performed  for  each 
value.  Each  finite  element  calculation  output  a  value  for  F  and  a  value  for  J  (energy  domain 
integral  method),  g  was  plotted  versus  a/w .  A  second  order  polynomial  was  fitted  to  g , 


=  P0  +  Pl^  +  P2 


(9) 


giving  a  close  approximation  to  the  function.  From  Equation  8, 
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Kl(\a,b,w,h)  = 

Using  Equations  7  and  9  in  Equation  10  gives 

K1(A,a,b,w,h  =  5)  =  J§"{p^  +  p^  +  P2{w))- 

For  each  experimental  point  (A,  a)  of  each  specimen  thickness  ( b ),  KX(A,  a,b,w  =  1)  was 
determined  using  Equation  11,  in  which  F  was  the  experimentally  measured  specimen  end  force 
associated  with  (A,  a)  . 

The  alternative  “brute  force”  method  for  calculating  Kx  is  to  make  one  finite  element  calculation 

for  each  experimental  point  (A,  a)  of  each  specimen  thickness  (b) ,  in  order  to  find  J  for  the 
point,  and  then  to  compute  Kx  with  Equation  5,  using  E  =  E0.  The  advantage  of  using  the 

method  of  this  section  is  that  once  the  approximating  function  of  Equation  9  is  determined,  Kx 
can  be  found  for  any  number  of  points  (A,  a)  without  doing  more  finite  element  calculations. 
Each  finite  element  calculation  of  J  is  not  trivial  with  respect  to  the  number  of  elements  required 
to  obtain  a  reasonably  accurate  result.  When  determining  the  function  g  of  Equation  8,  only 
enough  points  were  (finite  element)  calculated  to  provide  the  function’s  shape. 

This  paragraph  is  a  clarification  of  the  previous  statement  that  g  depends  only  on  the  two 
parameters  a/w  and  h/w .  The  experimental  SENT  specimens  were  symmetric  with  respect  to 
geometry  and  loading,  about  the  crack.  So  for  the  finite  element  calculations,  one  half  of  a 
specimen  was  modeled,  cut  thru  the  crack  centerline.  The  mesh  boundary  along  the  crack 
centerline  was  referred  to  as  the  bottom.  The  points  of  the  mesh  bottom  boundary  which  were 
ahead  of  the  crack  tip  were  constrained  from  movement  in  the  y  (perpendicular  to  the  crack) 
direction.  The  displacement  A  was  applied  in  the  y  direction  to  each  point  on  the  mesh  top 
boundary.  The  ends  of  the  experimental  specimens  were  glued  to  aluminum  end  tabs.  In  order  to 
simulate  the  glued  condition,  each  point  of  the  mesh  top  boundary  was  constrained  from 
movement  in  the  jc  direction.  For  this  case,  in  which  the  simulated  specimen  is  prevented  from 
free  lateral  contraction,  the  induced  stress  distribution  in  the  interior  and  on  the  boundaries  of  the 
specimen,  is  some  unknown  nonlinear  function  of  position  ( x ,  y)  and  of  Poisson’s  ratio  (v),  i.e., 
a  =  c(x,  y,  v) .  This  results  in  g  at  any  point  being  dependent  on  the  value  of  v .  (This  effect  of 
v  on  g  also  exists  if  load  is  applied  to  the  specimen  by  prescribing  tractions.)  Since  g  depends  on 
v ,  g  is  also  dependent  on  whether  plane  stress  or  plane  strain  is  used.  The  influence  of  the  value 
of  v  on  the  value  of  g  increases  primarily  with  decreasing  h/w ,  and  secondarily  with  decreasing 
a/w .  As  h/w  increases,  g  becomes  independent  of  the  value  v ,  and  also  of  the  type  of  loading 
(prescribed  uniform  edge  displacement  or  prescribed  uniform  edge  stress).  The  finite  element 
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calculations  of  g  which  were  made  for  the  experimental  specimens  of  this  paper  can  be 
confidently  applied  only  for  the  particular  case  of  a  plane  stress  SENT  specimen  with  v  =  0.499 
and  h/w  =  5 ,  subjected  to  prescribed  uniform  edge  displacement.  There  is  no  way  to  convert 
these  values  of  g  to  a  different  case. 

The  Crack  Growth  Resistance  Curve 

Figure  11  is  a  plot  of  K1  vs.  the  change  in  the  crack  length  (5a) .  These  curves  are  called  crack 
growth  resistance  curves  (R-curves),  and  are  a  measure  of  a  material’s  resistance  to  crack  growth. 
The  average  value  of  at  crack  initiation  (Kic)  is  47.8  lb.  /  in.1’5.  The  maximum  value 

recorded  for  Kx  is  81.4  lb.  /  in.1'5 .  The  condition  Kl<Klc  corresponds  to  crack  tip  blunting,  as 
is  shown  in  Figure  1.  The  overall  trend  is  for  Klc  to  increase  with  increasing  b .  (This  trend 
provides  additional  evidence  that  the  thicker  specimens  were  not  in  a  state  of  plane  strain.  If  it  is 
assumed  that  as  a  specimen’s  thickness  (b)  is  increased,  the  state  of  the  crack  tip  region 
transitions  from  plane  stress  to  plane  strain,  then  Klc  should  decrease  with  increasing  b ,  because 
a  material’s  crack  growth  resistance  is  smaller  in  plane  strain  than  in  plane  stress.) 

Analysis  of  the  Crack  Growth  Rate 

Figure  12  shows  a  plot  of  a  vs.  K1 .  The  distribution  of  the  points  suggest  the  use  of  a  parabolic 

modeling  function.  According  to  Knauss  (6),  and  Schapery  {7},  the  crack  growth  rate  for  a 
SENT  specimen  of  a  linearly  viscoelastic  material  can  be  expressed  as 

d  =  cxK\2.  (12) 

(Note  that  in  Equation  12,  Cj  and  c2  are  not  asserted  to  be  material  properties,  but  rather  are 
assumed  to  be  specimen  dependent.)  Therefore, 


log10(d)  =  logjotc^  +  c^ogjoCfiT!).  (13) 

Figure  13  shows  a  plot  of  log10(d)  vs.  log10(A’1) .  Also  shown  is  a  least  squares  linear  curve  fit, 
giving  (for  [a]  =  in.  / min.,  [ATj  ]  =  lb.  /in.1’5 ,  (  [jc]  =  "units  of  x ") ) 

log10(d)  =  -4.58  +  2.26  log1()(isri).  (14) 

.  .  -4  58 

Therefore,  the  best  fit  pair  (cj,  c2)  for  representing  all  eight  SENT  specimens  is  ( 10  '  ,  2.26) . 
Equation  12  can  then  be  written 


d  =  2.63x10  5^2'26 


in./min. 


(15) 
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It  was  observed  that  cx  and  c2  are  not  independent  parameters.  Figure  14  shows  a  plot  of 
logjoCcj)  vs.  c2  (one  best  fit  pair  (cv  c2)  is  calculated  for  each  of  the  eight  SENT  specimens). 
Also  shown  is  a  least  squares  linear  curve  fit,  giving, 

logjoCcj)  =  -0.553  -  1.79c2.  (16) 

Substituting  Equation  16  into  Equation  12  eliminates  c1  and  gives 


a  =  0.280  •  (0.01 64KX)C2  in./min. 
c2  can  be  solved  for,  giving 


log10(a)  +0.554 
2  =  log^)-!^- 


(17) 


(18) 


Equation  18  indicates  that  for  a  particular  SENT  specimen,  any  experimental  point  (Kx,  a) 
should  give  an  estimate  for  the  parameter  c2  for  that  specimen. 

Substituting  Equation  16  into  Equation  13  and  rearranging  gives 

log10(d)  + 0.553  =  c2  •  (log10(AT1)  -  1.79) .  (19) 

Equation  19  has  the  form  of  a  line,  i.e., 

y-yx  =  m-  (x-xx).  (20) 

Therefore,  Equation  16  is  identically  satisfied  for  any  line  which  passes  thru  the  point 
(logjoCATj^)  =  1.79,  log10(a)  =  -0.553),  and  has  some  slope  c2.  This  means  any  line  passing 

thru  the  point  (1.79,  -0.553) .  If  Equation  16  was  exactly  true  (if  each  point  in  Figure  14  was 
exactly  on  the  line),  then  in  Figure  13  the  shown  line  (which  is  the  best  fit  for  the  aggregate  of  the 
specimens)  would  pass  thru  the  point  (1.79,  -0.553) .  Also  if  Equation  16  was  exactly  true,  then 
in  Figure  1 3  if  a  best  fit  line  was  plotted  for  each  specimen,  all  of  those  lines  would  pass  thru  the 
point  ( 1 .79,  -0.553 ) .  Therefore,  point  ( 1 .79,  -0.553 )  can  be  designated  the  “pivot  point”  for  line 
fits  of  specimen  experimental  data  of  log10(d)  vs.  log10(12'1) . 

Each  point  (log10(cj),  c2)  in  Figure  14  was  determined  by  calculating  a  best  fit  line  thru  the 
points  (log10(7T1),  log10(a))  for  the  particular  specimen.  Then  the  relation  of  Equation  16  was 
determined  by  calculating  a  best  fit  line  thru  all  of  the  points  (log10(c1),  c2).  The  relation  of 
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Equation  14  was  determined  by  calculating  a  best  fit  line  thru  the  points  (log10(^ i  )>  logio(^) )  f°r 
all  of  the  specimens  together.  From  Equation  14,  log^Cj)  =  -4.58.  If  c2  from  Equation  14 
(c2  =  2.26)  is  substituted  into  Equation  16,  then  log jq(c j  )  =  —4.60.  The  two  values  of 
log10(cj)  vary  by  0.4%.  Their  agreement  is  good  because  each  point  in  Figure  14  is  close  to  the 
fit  line. 

Equation  17  can  be  written  in  the  form 

a  =  c3  ■  {caKx)2  in./min.,  (21) 

with  c3  =  0.280  and  c4  =  0.0164.  c3  and  c4  are  independent  of  specimen  thickness  (b),  and 
are  constants  for  this  particular  composite  solid  propellant,  and  the  particular  values  for  the 
specimen  width  (vv  =  1  in.) ,  specimen  length  (h  =  5  in.) ,  notch  length  (or  initial  crack  length, 
a0  =  0.3  in. ),  and  applied  strain  rate  (e  =  0.04  in./in./min. ).  c2  is  a  variable,  dependent  on  the 
specific  randomness  of  a  particular  specimen.  Equation  21  is  only  approximately  true,  because 
Equation  16  is  only  approximately  true.  For  an  actual  specimen,  when  Kx  equals  Klc,  then  a 

equals  0.  In  order  for  Equation  21  to  approximate  this  condition,  c2  should  increase  as  Klc 
increases.  For  any  value  of  c2 ,  a  curve  having  the  form  of  Equation  21  will  pass  thru  the  common 

intersection  point  ^Kx  =  a  =  c^j.  The  common  intersection  point  for  best  fit  (to  determine 

c2)  curves  of  each  of  the  eight  SENT  specimens  represented  in  Figure  12  would  be 

(Kx  =  61.0  lb.  +  in.1'5,  a  =  0.280  in./min. ).  At  the  common  intersection  point,  the  slope  of  the 
function  of  Equation  21  is  linearly  proportional  to  c2,  and  therefore  increases  with  Klc  (because 
c2  increases  with  Klc).  This  means  that  at  the  common  intersection  point,  the  rate  of  change  of  a 
with  respect  to  Kx ,  increases  as  KXc  increases. 

Conclusions 

Constant  uniaxial  strain  rate  tests  were  performed  on  four  thicknesses  (0.2  in.,  0.5  in.,  1.0  in.,  1.5 
in.)  of  specimens  of  a  particular  composite  solid  propellant.  Tests  were  performed  on  eight 
smooth  (unnotched)  specimens  to  determine  the  material’s  uniaxial  stress-strain  behavior.  Tests 
were  performed  on  eight  single  edge-notched  tension  (SENT)  specimens  to  determine  cracked 
uniaxial  stress-strain  behavior,  the  material's  Mode  1  crack  growth  resistance,  and  to  determine 
and  model  the  crack  growth  rate  (a).  Methods  were  described  for  calculating  a  from  crack 
growth  data,  and  for  calculating  the  Mode  1  stress  intensity  factor  (Kx).  For  all  analysis,  linear 
elastic  fracture  mechanics  (LEFM)  was  used.  Results  showed  that  there  were  no  consistent  and 
significant  effects  of  specimen  thickness  (b) .  The  crack  growth  data  made  it  seem  likely  that  for 
all  four  values  of  b ,  the  crack  tip  region  experienced  a  state  of  plane  stress.  This  was  explained  as 
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being  due  to  void  formations  in  the  crack  tip  damage  zone  preventing  the  lateral  constraint  which 
would  otherwise  produce  a  state  of  plane  strain  in  a  thicker  specimen.  The  crack  growth  rate  was 

modeled  as  a  =  .  cx  and  c2  were  assumed  to  be  mutually  independent  of  each  other,  and 

to  individually  depend  on  the  specific  randomness  of  a  particular  specimen.  Then  data  showed 
that  Cj  and  c2  were  mutually  dependent,  and  one  could  be  eliminated.  The  crack  growth  rate 

could  then  be  modeled  as  a  =  c3  ■  {cAKx  f2 .  c3  and  c4  were  independent  of  b ,  and  depended  on 
the  particular  material,  and  the  particular  values  for  the  specimen  width  (w  =  1  in.),  specimen 
length  (h  =  5  in.),  notch  length  (a0  =  0.3  in.),  and  applied  strain  rate  (e  =  0.04  in./in./min. ). 
So  c3  and  c4  were  constants  for  the  eight  SENT  specimens.  c2  depended  on  the  specific 
randomness  of  a  particular  specimen.  It  was  noted  that  all  SENT  specimen  curve  fits  (to 

determine  c2)  of  a  =  c3  •  (c^Kf)  2  would  pass  thru  the  common  point  =  l/c4,  a  =  c3), 
and  at  that  point  the  rate  of  change  of  a  with  respect  to  Kx  would  increase  with  increasing  values 

of*lc- 
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Figure  1 .  b  =  1 .0  in.,  Crack  Tip  Blunting  Figure  2.  b  =  1 .0  in.,  Initial  Crack  Growth 


Figure  3.  b  =  1.5  in.,  Damage  Formation  Figure  4.  b  =  0.2  in.,  Damage  Zone  Voids 


Figure  5.  b  =  0.2  in.,  Single  Ligament  Figure  6.  b  =  0.5  in.,  Double  Ligament 
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Figure  11.  K ,  vs.  8a 
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Figure  12.  a  vs.  Kx 


fit:  log10(a)  =  -4.58  +  2.261og10(AT1) 


log10(c,) 


fit:  log10(c,)  = -0.553 -1.79 
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